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ABSTRACT 

We present new, high-angular resolution interferometric observations with the Karl 
G. Jansky Very Large Array of 12 CO J=l — 0 line emission and 4-8 GHz continuum 
emission in the strongly lensed, z = 2.3 submillimetre galaxy, SMM J21352-0102. 
Using these data, we identify and probe the conditions in ~ 100 pc clumps within this 
galaxy, which we consider to be potential giant molecular cloud complexes, containing 
up to half of the total molecular gas in this system. In combination with far-infrared 
and submillimetre data, we investigate the far-infrared/radio correlation, measuring 
qm = 2.39 ±0.17 across SMM J21352. We search for variations in the properties of the 
interstellar medium throughout the galaxy by measuring the spatially-resolved qir and 
radio spectral index, a ra di 0 , finding ranges <?ir = [2.1, 2.6] and a ra dio = [—1.5,—0.7]. 
We argue that these ranges in a ra dio and qm may reflect variations in the age of the ISM 
material. Using multi-J 12 CO data, we quantitatively test a recent theoretical model 
relating the star-formation rate surface density to the excitation of 12 CO, finding good 
agreement between the model and the data. Lastly, we study the Schmidt-Kennicutt 
relation, both integrated across the system and within the individual clumps. We 
find small offsets between SMM J21352 and its clumps relative to other star-forming 
galaxy populations on the Schmidt-Kennicutt plot - such offsets have previously been 
interpreted as evidence for a bi-modal star-formation law, but we argue that they can 
be equally well explained as arising due to a combination of observational uncertainties 
and systematic biases in the choice of model used to interpret the data. 

Key words: galaxies: active — galaxies: high-redshift — galaxies: starburst — sub- 
millimetre — ISM: molecules — galaxies: ISM 


1 INTRODUCTION 

The resolution of the cosmic submillimetre background 
into discrete sources has revealed a population of high- 
redshift, ultra-luminous (Lm b 10 1 1 ~ 13 Lq) galaxies (e.g . 
Ilvison et al .1 1 1 998l : I Walter et al.ll2012l : iRiechers et al .1 [FoTtji . 
the submillimetre galaxies (SMGs), whose bolometric lumi¬ 
nosity is dominated by the re-radiance at infrared wave¬ 
lengths of dust-obscured UV/optical photons produced by 
young stars and, in some cases, from accretion on to a central 
super-massive black hole, seen as an active galactic nucleus 
(AGN). 


* E-mail: alasdair.thomson@durham.ac.uk 


In spite of the intrinsically high luminosities of SMGs, 
many remain challenging targets for observational studies, 
due both to photon starvation in the 1 mm atmospheric win¬ 
dow commonly used to survey for them and the low angular 
resolution of these single-dish surveys. Indeed, with the ex¬ 
ception of a handful of hyper-luminous objects (HyLIRGs), 
the brig htest SMGs have 850-um flux densit i es of only ~ 
lOmJy (Ilvison et al.ll2010bl : lKarim et al.ll2013l : lBarger et al] 
l2014li . These effects conspire to suppress many examples of 
this important class of star-forming galaxy at or below the 
detection limit of confusion-noise limited single-dish sur¬ 
veys, hindering efforts to understand how and where they 
form their stars within their interstellar medium (ISM). 
To alleviate these issues, submillimetre surveys have been 
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devised which exploit the panchromatic flux-boosting ef¬ 
fects of gravitational lensing, either by individua l li ne-of- 
sight foreground gal axies (e.g. iRowan-Robinson et ahlflQQll : 
iNegrello et alJ [2010) or by deliberate or serendipitous ex¬ 
ploitation of massive foreground clusters (ISmail et al.ll 19971 : 
iKnudsen et al.ll2008l : iKneibllioiCh . 

Among the most striking star-forming galaxies detected 
in this manner is “The Cosmic Eyelash” (SMM J21352-0102; 
hereafter SMMJ21352 at RA: 213511 . 6, De c: -0102 52, 
J2000), discovered by iSwinbank et~ahl (l2O10Tl following a 
survey of the z = 0.325 MACSJ2135-01 cluster with the 
Large Apex Bolometer Camera (LABOCA, ISiringo et ali 
l2009h on the 12-m Atacama Path-finder Experiment 
(APEX) telescope. The high observed flux density of 
SMM J21352 - SWo |im = (106 ±7) mjy - is a result of the ~ 
37.5 x magnification provided by the foreground cluster, and 
corresponds to a de-boosted flux density that is close to the 
confusion limit of submillimetre surveys (SWo^m ~ 3mJy). 
Its redshift, z = 2.3259, derived through the first blind de¬ 
tection of 12 CO J=l —0 with the Zpectrometer instrument 
on the Robert C. Byrd Green Bank Telescope (GBT) is 
close to the median redshift (z ~ 2.5) of the SMC population 
IChapman et alfeoosl : ISimpson et al.ll2014h , suggesting that 
SMMJ21352 may be an intrinsically typical SMC, located 
serendipitously behind a strong cosmic lens. In addition to 
boosting its total flux above the confusion noise, the high 
magnification of SMM J21352 spatially extends the emission 
in the image plane, allowing individual star-forming regions 
within this high-redshift starburst galaxy to be identified 
and studied in detail. 

Following the detection of SMMJ21352, 
iSwinbank et al] (|201C)h conducted a follow-up programme 
of observations with the Submillimeter Array (SMA) at 
870 /im, where the 0.3” synthesized beam in the very 
extended (VEX) configuration resolved the emission in to 
eight distinct clumps over a 4” region, representing four 
unique dusty clumps of ~ 100 pc size each (in the source 
plane), mirrored about the critical curve. 

The bulk molecular gas properties of SMM J21352 were 
first studied in [Danielson et al.l (l201lh , in which a combi¬ 
nation of single-dish 12 CO J = 1 — 0, mid-J 12 CO, Cl and 
HCN observations from the GBT, IRAM 30-m telescope 
and Plateau de Bure Interferometer (PdBI) were used to 
constrain the source-averaged excitation in SMMJ21352. 
These authors have found evidence for a two-phase medium 
in SMMJ21352 comprising a hot, dense, luminous compo¬ 
nent embedded within a mass ive (> IQ 10 Mp>) , exten ded, 
cool, low-excitation reservoir. iDanielson et al.1 (j201lh use 
their spatially-integrated 12 CO spectral line energy distri¬ 
bution (SLED) and spectrum to fit a kinematic model com¬ 
prising four distinc t clumps, with d if fering physical proper¬ 
ties. Subsequently, ISwinbank et all ( 2011 ) studied the gas 
kinematics using moderate-resolution (~ 0.8”) observations 
from the Karl G. Jansky Very Large Array (VLA) of 12 CO 
J=l —0 emission, determining that they are well-described 
by a rotationally-supported disc model with an inclination- 
corrected rotation speed u ro t = 320±25 kms -1 , a rotational- 
to-dispersion velocity ratio v/a = 3.5 ± 0.2 and a dynam¬ 
ical mass of (6.0 ± 0.5) x IQ 10 Mp within a 2.5 k pc radius 
in the source plane. Later, iDanielson et al.l ( 2013 1 revisited 
SMM J21352 in a study of the optically-thin 13 CO and C 18 0 
lines, finding evidence of a C ls O over-abundance (relative 


to the Milky Way) in high-temperature (> 100 k), dense 
(n > 10 3 cm -3 ) regions, which may favour a somewhat top- 
heavy initial mass function. 

Here, we present the latest high-resolution, high- 
sensitivity interferometric 12 CO J = 1 —0 map and spec¬ 
tral data cube of SMM J21352 obtained with the upgraded 
VLA, together with a new, deep VLA C-band continuum 
and spectral index map. We first use these observations 
to spatially-identify, and then study the excitation condi¬ 
tions within dense molecular gas clumps, providing an im¬ 
portant independent check on previous results derived via 
kinematic decomposition of the source. In these clumps, we 
carry out the first observational test of recent theoretical 
models which relate the 12 CO SLED directly to the star for¬ 
mation rate (SFR) density Esfr- Next, we use our new data 
to create a spatially-resolved Schmidt-Kennicutt (SK) plot, 
to search for evidence that either the clumps in SMM J21352 
or the bulk ISM deviate strongly from the star-formation law 
ascribed to normal galaxies. We also study the spatially- 
resolved far-infrared/radio correlation (FIRRC), to search 
for evidence of possible age variations within the ISM. 

In § 2, we present our observations and discuss the 
data reduction strategies for our new VLA data. In § 3 
we present the latest constraints these data provide on the 
source-averaged properties of SMMJ21352. §4 utilizes the 
high angular resolution of our observations to measure the 
gas excitation conditions, radio properties and ratio of far- 
infrared/radio emission within the ISM, and in § 5 we dis¬ 
cuss the implications of these results for our understanding 
of star-formation at high redshift, and give our conclusions. 

Throughout this paper, we use a cosmology with Ho = 
71 kms -1 Mpc -1 , = 0.27 and Da = 0.73, which gives 

an angular scale of ~ 8.32/37.5 = 0.22 kpc arcsec -1 . All flux 
densities quoted are observed, and all luminosities - unless 
otherwise stated - are corrected for the 37.5x amplification. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 VLA 12 CO J=l — 0 observations 

Observations of redshifted 12 CO J = 1 — 0 emission in 
SMM J21352 were carried out using the VLA’s new A'a-band 
receivers in 2-4 hr blocks (totalling 25 h of usable on-source 
data) in excellent w eather conditions betwe en 2010 Septem¬ 
ber - 2011 May. In ISwinbank et ahl (|201lh we presented a 
datacube and map made from the hybrid DnC/CnB ob¬ 
servations that were available at that time (VLA projects 
10A-224 and 10B-195), in which the 12 CO J = 1 — 0 emis¬ 
sion was shown to be resolved on ~ 0.8” scales - here we 
incorporate the latest BnA-configuration (10C-141) obser¬ 
vations, in order to maximize our angular resolution and 
study the ISM on hitherto inaccessible scales. During this 
Open Shared Risk Observing (OSRO) period, the Wide¬ 
band Interferometric Digital ARchitecture (WIDAR) corre¬ 
lator provided 64 x 2 MHz dual-polarization channels in each 
of two sub-band pairs, corresponding to a velocity resolution 
of ~ 15kms -1 channel -1 before smoothing. 

Antenna pointings were checked hourly using the C- 
band receivers, and immediately prior to scans of the pri¬ 
mary flux calibrator, 3C48. Amplitude and phase corrections 
as well as the bandpass were determined via short scans of 


© 2014 RAS, MNRAS 000.IHIT51 
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Figure 1. (a): Velocity-inte grated VL A 12 CO J = l —0 map (red contours), showing the position of the critical curve from the lensing 
model (blue contour; ISwinbank et al . 201C ). The underlying greyscale is from the If ST /WFC3 F140W observations of the foreground 
cluster, MACS J2135-0102 iSwinbankeT^LbOllh . Contours are spaced in intervals of a = 14/My, beginning at 3<r. We mark with yellow 
dots the positions of the gas clumps identified in the 12 CO J = l —0 datacube via two independent m ethods f S I4.lt . and identify them 
as spatially-resolved counterparts to the kinematic clumps X, Y, Z\ and Z 2 of [Danielson et ahl j20 1 il l (b): Greyscale showing spatial 
variations in the C-band radio spectral index (o ra dio)- The red contours represent a spatially-smoothed version of the 4-8 GHz data 
presented in the fourth panel, (c): Updated SMA 870 (im image, made usin g prev iously existing uv data, but incorporating additional 
short spacings with respect to the VEX map presented in ISwinbank et al.l j2010l ). Contours are plotted in intervals of a = 1.4 mjy, 
beginning at 3<r. (d): High-resolution VLA C-band (4-8 GHz) continuum image, with contours starting at 3cr (where a = 1.8/tJy) and 
in steps of v/2 X a thereafter. In all four maps, the beam used to create the contour image is shown at the bottom right. 


the bright, nearby calibration source J2136+0041 in between 
each 5 min integration on SMM J21352. 

The data were flagged, calibrated a nd imaged using th e 
standard ALT'S procedures outlined in llvison et alj (20 11). 
though with a number of important changes: the data were 
downloaded in native SMBDF format and loaded into ALT’S 
using bdf2aips rather than as uvfits files using fitld, in 
order to eliminate any server-side data compression arte¬ 
facts. Known issues in the setting of delays at the VLA were 
compensated for in each track using fring, and 1 min of 
data towards the end of a good scan of 3C48. 

For the DnC observations, the two sub-bands were 
tuned to 34.62 and 34.73 GHz, such that at the GBT red- 
shift (2 = 2.3259; ISwinbank et aLl |2010l) t he 12 CO -J = 1—0 
line (u rest = 115.27120256 GHz; iMorton fc Noreaul Il994l ) 


ought to have been centred in channel 50 of the lower sub¬ 
band (AC), with the upper sub-band (BD) tuned to overlap 
by 10-channels to mitigate the effects of noise in the end 
channels of each sub-band. The resulting 118 channels pro¬ 
vide 2050 km s -1 of usable velocity coverage. In fact, spectra 
from these DnC-only data hinted at a slight redshift discrep¬ 
ancy with respect to the reference GBT redshift, and hence 
the array tunings for the subsequent CnB and BnA obser¬ 
vations were shifted downward in frequency by 40 MHz to 
accommodate the full width of the line. Thus, the usable 
spectral coverage common to all three array configurations 
is 98 channels, covering 1690 km s _1 . 

Calibrated uv data for the three array configura¬ 
tions were concatenated together with dbcon, with the 
reweight parameter set to provide a 1:4 ratio in the sum 
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of gridded weights between adjacent configurations, and so 
maximize the likelihood of the final data set producing a 
Gaussian synthesized beam when imaged. 

We created maps and spectral data cubes of the 12 CO 
J = 1 — 0 emission over a 3-arcmin field of view, which we 
cautiously cleaned by placing small, circular clean boxes 
around the brightest visible knots of emission after each 
major CLEAN cycle, then inspecting the residual image and 
iteratively adding more small clean boxes as appropriate un¬ 
til the network of small clean boxes overlapped across the 
full extent of the source and all emission lying significantly 
above the noise had been cleaned. Such conservative cleaning 
is necessary when creating interferometric maps of low sur¬ 
face brightness emission in order to prevent spurious emis¬ 
sion (e.g. noise, or sidelobe emission) being carried through 
to the final CLEAN image, where it will lead morphological 
analysis astray. 

We applied Gaussian uv tapers from 50-1000 kA to the 
data in steps of 50 kA to search for any flux which may have 
been resolved out on long baselines, but found no evidence 
of significant missing flux. The un-tapered map (Figure, [l]) 
contains 2.07 ±0.01 mjy, ~ 95% of the single-dish GBT flux 
llSwinbank et al.ll2010ll . and hence we use this map without 
any additional tapering. The synthesized beam is 0.28 x 
0.21 arcsec at position angle PA= —50°. 

2.2 VLA C-band continuum observations 

SMM J21352 was observed in the extended A-configuration, 
in C-band (4-8 GHz) at the VLA during 2011 July 6 th -7 th , 
in moderate weather conditions. These observations took 
place towards the end of the (Expanded) VLA hardware 
upgrade phase, by which time all 27 antennas had been fit¬ 
ted with functioning C-band receivers, and the WIDAR cor¬ 
relator was tuned to provide 16 sub-bands of 64 x 2 MHz 
channels each, yielding a total instantaneous bandwidth of 
2 GHz. By splitting the allocated time between observations 
at 4-6 and 6-8 GHz, a moderately large spectral region of 
the radio continuum emission in SMM J21352 was studied. 

From a total time allocation of 10 h, 8 h of usable on- 
source data were obtained, with flux calibration being set via 
scans of 3C48 at the end of each track, and antenna pointing, 
amplitude and phase tracking and bandpass correction being 
determined via frequent, short scans of the nearby calibrator 
J2136±0041. 

These data were reduced and imaged in the Common 
Astronomy Software Applications (CAS40) package. Using 
the MSMFS CLEAN algorithm, we produce a high-resolution 
map with a synthesized beam 0.34” x 0.29” at PA= 180° 
and an rms sensitivity of 1.45 pjy beam -1 . The total flux 
density measured in this map is Sc = 0.258 ± 0.001 mjy. 
This map is shown in FigureQ] 

2.3 IRAM PdBI mid-J 12 CO observations 

We use observations from the six-element IRAM PdBI to 
study the mid-J 12 C O emission. These obse rvations were 
originally presented in lDanielson et ahl (120111 1. Using 4 and 

1 http://casa.nrao.edu 


2 hr on-source exposures in D-configuration, taken in Di¬ 
rectors Discretionary Time (DDT) during 2009 May, 12 CO 
J = 3 — 2 and J = 4 — 3 emission were observed, reach¬ 
ing signal-to-noise (S/N) ~ 300 across 900kms _1 in each 
velocity-integrated map. The total measured line fluxes are 
13.20 ± 0.10 Jy kms -1 and 17.3 ± 1.2 Jy km s -1 for J = 3 — 2 
and J = 4—3 respectively, and the corresponding synthesized 
beams are 2.11” x 1.59” at PA= 37° and 1.08” x 0.50” at 
PA= 23°. Furt her details of the data reduction and imaging 
are presented iDanielson et ahl (j201lf h 

12 CO J = 6 —5 emission was observe d in A c onfigura- 
tion during 2010 January fsee lSwinbank et al.ll201 ill . The six 
hours of usable on-source integration time in this configura¬ 
tion yield a synthesized beam of 0.67” x 0.43” at PA= 24.5° 
and a velocity-integrated line flux of 21.5 ± 1.1 Jykms -1 . 

2.4 SMA observations 

A programme to map the rest-frame far-infrared continuum 
emission in SMMJ21352 w as carried out at the SMA in 
2009 May - September (see ISwinbank et al.ll20ld . in which 
a description of the data reduction is presented). With the 
lower sub-band tuned to 870 /im, a total of 5.5 h of on-source 
data were obtained in the sub-compact (SC) configuration, 
with 12 h in compact (C), 2h in extended (Ext) and 6h in 
the very extended (VEX) configuration. The upper sub-band 
was tuned to search for redshifted 12 CO J=10 —9 emission 

- this line was not detected, but a formal 3<r upper-limit of 
< 0.2 Jy kms -1 was determined. 

In order to search for signs of clumpy dust emission, 
ISwinbank et alJ (hoiol ) used only the six hours of VEX- 
configuration uv data to produce a high-resolution 870 /tm 
image. The resulting map - with a 0.3” synthesized beam 

- comprises eight, bright image-plane clumps (four source- 
plane clumps, each mirrored about the critical curve of the 
foreground lens), which account for ~ 85% of the total (i.e. 
single-dish LABOCA) 870 pm flux. In light of the updated 
information on the morphology of this system due to the new 
VLA A'a-band (covering 12 CO J = 1—0) and C-band (contin¬ 
uum) imaging, we opt to carry out a new imaging run from 
the existing calibrated SMA 870 pm visibilities. We find that 
by incorporating previously unused visibilities from the Ext, 
C and SC array configurations - though with a weighting 
scheme which still favours the longest baselines - we suffer 
some degradation in resolution (0.59” x 0.45” at PA = 50°, 
cf. 0.29” x 0.23” at PA = —15° from VEX-only), but with 
the advantage that the improved sampling of the uv plane 
yields a synthesized beam with less intrusive sidelobes. This 
new, more conservative SMA image (Figure 1) has a sensi¬ 
tivity of a = 1.4 mjy beam -1 , and the total flux density of 
5*870 /im = 91 ± 5mJy acco unts for ~ 85% of the single-dish 
LABOCA flux reported in lSwinbank et ahl d2010l l. 

3 SOURCE-AVERAGED PROPERTIES 

3.1 12 CO J=l —0 spectral line profile 

SMM J21352 is detected in 12 CO J=l —0 with a peak flux 
in the velocity-integrated image of 8cr, where a = 14pJy. 
We create a spectral data cube with 6 MHz spectral res¬ 
olution by binning the native VLA channels in groups of 
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Figure 2. SMM J21352 12 CO J = 1—0 and J — 3—2 spectra. The 
latter is scaled vertically by a factor 9 —1 so as to lie on the same 
Tb scale. The velocity scale is quoted with respect to the GBT 
12 CO J = 1 — 0 redshift, and a typical 50//Jy/channel error bar 
is shown in the top-right of the plot. The shapes of the two line 
profiles are in good agreement. The new VLA J = 1—0 line profile 
recovers ~ 95% of the flux of t he GBT Zpectrometer J = 1 — 0 
measurement (not shown) from ISwinbank et alJ <|2010h . but has 
a slightly different line profile due to the superior velocity resolu¬ 
tion of the VLA’s new Ka -band receivers relative to the coarser 
channels of the much-wider bandwidth Zpectrometer instrument 
on GBT. We fit this new 12 CO J= 1 — 0 spectrum with emission 
from four kinematic components, which are well-modelled as a 
series of Gaussian components ( X , Y and Z± and Z 2 : faint grey 
dotted lines). 


three in order to maximize signal-to-noise, and extract a line 
spectrum across the source using the AXVS task blsum. 
We see evidence of a broad line in the spectrum (fwzi= 
920 k ms" 1 ; in ag ree ment w ith the existing mid-J 12 CO data 
from [Danielson et al Eon]), with the velocity-integrated flux 
being consistent with that measured in the equivalent inte¬ 
grated image. This line profile places SMM J21352 at a flux- 
weighted redshift of z = 2.3268 ± 0.0001. The uncertainty 
in this redshift corresponds to velocity uncertainties of the 
order AV~ 10kms -1 . We produce an error spectrum using 
the RSPEC task, in an off-source region of the cube. The typ¬ 
ical noise in each 6 MHz spectral channel is 51 pJybeam -1 . 

The peak J = 1 — 0 flux density seen in the source- 
averaged spectrum (Figure[5J) of ~ 7 mjy is slightly higher 
th an the peak in the sin gle-dish GBT spectrum presented 
in ISwinbank et al.l (l20ld f , yet the velocity integrated flux 
densities agree within 5%. We attribute the differences in 
line shapes between these line profiles to the ~ 3x higher 
velocity resolution of the VLA data cube, and note that the 
new J = 1 — 0 line profile is similar t o the shape of th e 12 CO 
J = 3 — 2 line profile presented in [Danielson et al . I (l201l|, 
,which we show in Figure[2]for comparison), in which it is 
found that emission is better-fit by a superposition of Gaus¬ 
sian components with velocity offsets relative to the system¬ 
atic zero velocity for SMM J21352 than they are by a single 
Gaussian component. Motivated by this analysis, we fit four 
Gaussians - representing a possible interacting system - over 
the fwzi= 920km s -1 width of the 12 CO J =1 — 0 line. The 
peaks and central velocities of these four components are 
given in Tablc[l] 


Table 1. SMM J21352 kinematic parameters 


Component 

Peak 

(mjy) 

V 

(kms -1 ) 

X 

2.4 

+202 ± 30 

Y 

6.7 

+ 120 + 30 

Zx 

2.0 

-70 ± 30 

Z 2 

1.8 

-310 + 30 


Notes: Best-fitting coefficients for the four-Gaussian fit to the 
observed 12 CO J = 1 — 0 line. Velocities are quoted with respect 
to the reference redshift of z = 2.3259, and uncertainties in the 
velocities are calculated by adding the statistical errors of the fit 
and the half-channel width in quadrature. 

3.2 Line luminosity and total gas mass 

Following lSolomon et al.l ( 1997 1. the intrinsic (i.e. de-lensed) 
12 CO J =1 — 0 line luminosity is given by 


L'r 


3.25 x 10 7 


CO(l-O) 


D 2 (Mpc) 


Jav S ’' dv 

Jy kms -1 


1 + Z 


K km s 1 pc 2 


^rest 

GHz 


(i) 


where fi is the average magnification of the source, Dl is 
the luminosity distance (Mpc), i+est is the rest-frame cen¬ 
tral frequency of the line and f AV Svdv is the velocity- 
integrated line intensity. We measure Lco(i-o) = (1.7 ± 
0.2) x 10 10 Kkms -1 pc 2 - in the limit that the molecu¬ 
lar gas is optically thin, this implies a lower-limit on the 
(galaxy-averaged) 12 CO-to-H 2 conversion factoiQ_aco_> 
0.54, and hence M gas > 9.1 x 10 9 Mg (see llvison et al.1 
1201 ll : iDanielson et_ak[ 2011]^ 2013|) , whi le adopting the dy- 
namical model of ISwinbank et all (1201 il l and assuming that 
the dynamical mass is entirely dominated by molecular gat(3 
gives corresponding upper-limits aco < 3.5 and M gas < 
6.5 x 10 10 M b . 


3.3 Radio spectral index 

As a result of our strategy of splitting our VLA C-band 
observations between two different central frequencies, our 
observations cover a total bandwidth A v = 4 GHz, com¬ 
prising 32 IFs of 64 x 2 MHz channels. This relatively large 
fractional bandwidth provides a long enough lever along 
which to measure the local radio spectral index. Setting the 
CASA CLEAN parameter nterms= 2 implements the multi¬ 
scale, multi-frequency synthesis (msmfs) algorithm, which 
uses two Taylor terms to model the sky emission, and si¬ 
multaneously produces a cleaned continuum image, spectral 


2 «co includes a ~ 36% correction, based on primordial nucle¬ 
osynthesis, for the contribution to the ISM of the second most 
abundant element, He. The molecular hydrogen mass is therefore 
Mh 2 ~ Mg as /1.36. 

3 While such a scenario is physically implausible, we note 
the difficulty in accounting for the stellar m ass in SMGs (see 
IWardlow et al.l|201ll : lMichalowski et a.l~l1 201 T) . and consider that 
our “gas dominated dynamics” toy model at least gives a firm 
upper limit on M gas . 
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index (a ra dio) image, and spectral index error (da ra dio) im¬ 
age. The uncertainty, da ra dio as a function of position in 
SMM J21352 is inversely proportional to the signal-to-noise 
in the continuum image - we mask our spectral index image 
by excluding all pixels with |da ra dio| > 0.2, which corre¬ 
sponds to a S/N cut of ~ llcr in the continuum image. 

The median spectral index of all pixels within this error 
threshold is a ra dio = —0.96 ± 0.13, which is consistent with 
th e integrate d Ka-, X - and C-band flux densities measured 
in llvison et al.l ll2010all . which imply a rad i 0 = —0.22 ± 0.63. 

The spectral index map is shown in Figure[T] with the 
positions of the 12 CO J = 1 — 0 clumps marked for refer¬ 
ence. We search for signs of variations in a ra dio with posi¬ 
tion along one dimension by collapsing the spectral index 
(and error) image(s) along an axis aligned with the posi¬ 
tion angle of SMM J21352 and measuring the median o? ra dio 
and da? ra dio in each resolution step along this axis. We mea¬ 
sure the spatially-reoslved gas velocity field by extracting 
a position-velocity (PV) diagram from the 12 CO J = 1 —0 
datacube along the same slice, using the CASA task impv. 
The Q ra dio profile as a function of position from the critical 
curve for the brighter SMM J21352 counter-image is shown 
in Figure[31 We see that radio emission located 0.7-1.5” from 
the critical curve is well described by a ra dio ~ —0.7 (consis¬ 
tent with non-thermal synchrotron emission associated with 
star-formation), but that closer to the critical curve the radio 
spectrum noticeably steepens (to a ra dio ~ —1.5). It is plau¬ 
sible that this region of steep-spectrum emission is the result 
of ongoing merger activity, where two (or more) star-forming 
progenitors spawn radio continuum bridges and tidal tails as 
the interaction proceeds. Such tidal tails are not associated 
directly with star-formation activity, but may produce ultra¬ 
steep spectrum (a ra dio < 1) synchrotron emission as a re¬ 
sult of relativistic electrons being shock-accelerated by the 
merger and confined within the in terpenetrating magnetic 
fields of the progenitor galaxies dLisenfeld &; Volkl l20ld : 
lMurphvll2013l l. We attempted to extract the radial a? ra di 0 
profile of the fainter counter-image as well, but found that 
the lower signal-to-noise of the data in this image rendered 
the results inconclusive except in the most strongly-lensed 
~ 0.5” near the critical curve (where the results were con¬ 
sistent with, albeit with larger uncertainty than the results 
shown in Fig. [3]). 


4 SPATIALLY RESOLVED PROPERTIES 

4.1 Identification of molecular gas clumps 

By fittin g a kinematic model to their source- integrated 12 CO 
spectr a, iDanielson et al. | teoiil, 2013) and ISwinbank et al.l 
l|2011f l discovered that much of the molecular gas in 
SMM J21352 resides in dense clumps in the galaxy. Initially, 
these were identified with clumpy s t ructu res seen in the 
SMA 870/an map of ISwinbank et ahl J20ld l. which appear 
to account for ~ 80% of the total single dish flux. With our 
latest high-resolution 12 CO J = 1 — 0 datacube, we check 
these positions using two independent methods. 

First, we use the ALT’S task SERCH, whi ch employs 
the matched filter analysis technique outlined bv lUson et al.l 
dl99lli . and attempts to fit Gaussian line profiles of user- 
specified channel width to each pixel in the datacube. Em- 
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Figure 3. A plot showing variations in the radio spectral index, 
Q r a dio as a function of distance ( R ) from the critical curve for 
the brighter of the two counter-images of SMM J21352. Shown 
at the top of the figure is a cut-out of the radio spectral index 
map for this counter-image, with greyscale running from black 
(a ra dio = —2.0) to white (a ra dio = 0.0). The solid black line and 
shaded grey polygon represent the median and Her scatter in 
o ra dio (respectively) as a function of distance from the critical 
curve. The dashed line represents the velocity field (V), derived 
from the 12 CO J = 1—0 position-velocity diagram, while the dotted 
line represents (in arbitrary units) the variation in the strength 
of the magnifying effect with distance from the critical curve. 
We find that the bulk of the radio emission in SMM J21352 has 
“radio ~ —0.8, consistent with a star formation origin, but that a 
highly-magnified pocket of emission ~ 0.5” from the critical curve 
has considerably steeper spectral index, which may be the result 
of the on-going coalescence of two or more distinct components. 



ploying a 4<r S/N cut and searching for lines between 50- 
950 kms -1 (in steps of 100kms” 1 ), we identify eight bright 
clump features in the image plane (four to the north of the 
lens critical curve, and four mirror-images to the south), 
and use constraints provided by the lens model and by 
the velocity field to pair the eight image-plane clump fea¬ 
tures to four source-plane gas clumps. These gas clumps 
have slight positional offsets relative to the SMA clumps of 
ISwinbank et al.l i 201 ( :;, but do have velocity centroids (rela¬ 
tive to the systemic velocity of SMM J21352) that alow them 
to be identified with t he kinematic clumps X, Y, Z\, Zi of 
llDanielson et al1l201ll l. 

The robustness of these clumps is verified by performing 
a blind s earc h for clum ps using the idl clump-finding algo¬ 
rithm of iWilliams et ahl (i 19941 n. which contours the data 
cube at 3<r (where a is the typical rms per spectral channel, 
~ 50 /ijy beam -1 ) to identify peaks, and then attempts to 
follow the surrounding emission through adjacent beams un¬ 
til its extended structure attenuates below the noise thresh¬ 
old or blends with neighbouring clumps. We impose no prior 

4 http://www.ifa.hawaii.edu/users/jpw/clumpfind.shtml 
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image-plane or kinematic constraints on the algorithm, and 
the fwhm of each clump line profile is computed on-the-fly, 
as the algorithm examines the cube. This algorithm reports 
18 possible clumps in total, the brightest four of which agree 
positionally to within ±0.5x the synthesized beam with the 
four northern clumps identified via SERCH. We identify four 
likely counterparts on the other side of the critical curve, 
within the same positional uncertainty. The positions of 
these clumps are shown in the main panel of Figure[T] 

Though the latest, high-sensitivity (but lower resolu¬ 
tion) SMA 870 fjan. map fails to recover the clumpy structure 
seen in the previous VEX SMA map, it is noteworthy that 
we do observe clumpy structure in the 12 CO J = 1 — 0 dat- 
acube w hose kinematic propert ies a re similar to tho s e dis- 
cussed in lSwinbank et all (l20ld f and lDanielson et al.l (1201 ll . 
whose analysis was guided by the VEX SMA image), albeit 
with positional offsets between the 12 CO and SMA-identified 
clumps of the order of ~ 1 x 0.3” beam. Also noteworthy is 
the fact that these newly-identified clumps contain ~ 10% 
of the 12 CO J = 1 — 0 line flux, with the majority of the 
cold gas residing in a more diffuse structure within which 
the clumps are embedded. 

If these 12 CO gas clumps are marginally resolved by 
the VLA synthesized beam, then their physical sizes (once 
corrected for lensing) are ~ 100 pc, making them potential 
analogues of giant molecular clouds (GMCs), the sites of 
star formation in local galaxies. Do the physical conditions 
within the clumps justify this interpretation? 

4.2 Molecular gas excitation on ~ 100 pc scales 

We extract spectra from the 12 CO J = 1 —0 data cube 
at each of the clump positions. The resulting four clump 
spectra (Figure 2]) show velocity offsets with respect to the 
systematic velocity of SMM J21352. These offsets comprise 
one clump with a systematically positive velocity (which 
w e identify with the kin ematically-resolved component X 
of lDanielson et al.ll201lh : the next clump, with a pseudo- 
Gaussian velocity profile with narrow fwhm and broad fwzi 
(which matches component Y in the same paper); the re¬ 
maining pair of clumps have systematically negative veloci¬ 
ties (corresponding to clumps Z\ and Z 2 ). 

The physical conditions of the gas within these clumps 
are probed by first determining the 12 CO SLED of each 
clump, and then modelling the observed line ratios to con¬ 
strain the temperature, density and molecular abundance in 
each clump. 

We begin by placing the 12 CO J = 6 — 5, J = 4 — 3, 
J = 3—2, and J=l—0 maps on a common astrometry/pixel 
scale using the AXVS task OHGEO, and deblending the flux 
in the mid-J maps using a modified version of the algorithm 
presented in lSwinbank et ahl (|2014| j. Briefly, we seed model 
maps for each 12 CO transition with eight delta functions at 
the locations of the clumps (and their mirror images) iden¬ 
tified via our clump searching algorithms. The four clumps 
in the brighter mirror image north-west of the critical curve 
are constrained to have random flux densities between 0- 
150% of the peak seen in the corresponding real integrated 
image, and the flux densities of the four fainter mirror im¬ 
ages are constrained relatively via the lensing amplification 
map. Then, each model map is convolved with the appropri¬ 
ate synthesized beam, and the goodness-of-fit of that model 


map to the observed map is measured via the x 2 statistic, 
where the four degrees of freedom are those of the four in¬ 
dependent clump fluxes: 

2 _ Apix w X(Image - Model) 2 
X ~ 4 X E(Model) 2 1 ’ 

where N p i x is the number of pixels in each image (in each 
case, a 129 x 129 pix postage stamp with a pixel scale of 
0.04-arcsec). 

After 500 fits (one major cycle), we determine the set of 
clump fluxes for which y 2 reaches a minimum, re-centre on 
the best-fitting solutions and shrink the parameter space by 
5%, whereupon another set of 500 random clump fluxes are 
fit within the available parameter space. This process con¬ 
tinues for 40 major cycles (20,000 fits in all), or until Xmin 
has converged to within 1% across three successive major 
cyclefl 

Having thus determined the best-fitting 12 CO SLEDs 
for each of the clumps, these are then used to constrain a 
large velocity gradient (LVG) model, computed usin g the 
radex radiative transfer code (Ivan der Tak et ahl 120071 1. 

radex computes the level populations, optical depths 
and intensities of each of the spectral lines of a chosen set 
of molecules under the assumptions that: (i) the medium is 
homogeneous; (ii) that collisional excitation of the species 
occurs only via interactions with H 2 ; and (iii) that the ve¬ 
locity widths of all molecular lines are equal for a given 
cloud. The collision rates used by radex are functions of 
temperature, are held, along with other physical properties 
of each molecule, in a database of molecular hies maintained 
by the Leiden Atom ic and Molecular Database (LAMBDA0; 
ISchoier et al.1 120059 . and may be interpolated or extrapo¬ 
lated for Tk beyond the range contained within the datafile. 

A suite of 640,000 model 12 CO spectra were generated 
for each clump (with the line width altered appropriately 
in each case), for a grid of input parameters Tk=15-90 K 
(in 16 steps of 5 k), n(H2)=10 2-6 cm~ 3 (in 200 equidis¬ 
tant steps) and Nco = 4 x fo 18-24 cm -2 (in 200 equidis¬ 
tant steps), values typical of molecular clouds seen in the 
local Univer se (e. g. iRichardsonl Il985l : lOstriker et al.l 120101 : 
IShettv et al.ll201lh . In each model run, the background tem¬ 
perature was fixed at Tomb = 2.73 x (1 ± z) = 9.08 k, and 
the line-widths AV were set to the fwzi of the individual 
clumps (Figure[5]|. 

The measured 12 CO SLEDs for each of the four clumps 
are shown with the corresponding best-fitting radex mod¬ 
els in Figure[S] The shapes of the individual SLEDs are 
qualitati vely consistent with t he source-averaged SLED 
seen in 1 Danielson et al.1 d201lh . which may be a conse¬ 
quence of the assumption that the bulk of the excited 

5 In so doing, we implicitly assume that the emission from the 
./upper > 3 states is confined to the dense clumps - while this 
assumption is likely somewhat simplistic, we note that the high 
excitation conditions of the lines above J = 3 (n(H 2 ) ^ 10 4 cm -3 , 
Tj, y 50 k) suggest that they preferentially trace star-forming 
regions, and render it unlikely that a substantial contribution 
to this emission arises from the diffuse ISM. To overcome this 
assumption will require future high-resolution observations of the 
mid-J lines with ALMA. 

6 http://home.strw.leidenuniv.nl/~moldata/ 
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Figure 4. 12 CO J= 1—0 spectra for each of the four clumps identified via the methods presented in m The velocity offsets seen between 
ou r spatially-resolved clu mps allow them to be identified as likely counterparts to the kinematic clumps X , Y, Z\ and Z 2 , respectively, 
of I Danielson et all <|2Qlll ). The measured V'fwhm of each clump line profile is shown. 


gas phase is contained within these clumps. The SLEDs 
peak for J up per ^ 5, in keeping with results from stud- 
ies of local starburst galaxies (e.g. M82 and NGC6240; 
Weifi et al.||2005|: Meiier ink et al.ll2013l l. AGN (e.g. Mrk213; 


van der Werf et ah 2010T ). and from the few high redshift 


SMGs thus-far to have y ielded well-sampled SLEDs (e.g. 
GN 20 and SMMJ14011; ICarilli et alJ l2010l : ISharon et al.l 
120131 . respectively). The large error bars on the observed 
SLEDs are chiefly the result of degeneracies in the deblend¬ 
ing process, which arise due to the poor angular resolution 
of the mid- J PdBI maps. 


From the best-fitting radiative transfer models for each 
clump, we compute Xco = Nh2/7co, the ratio of molecu¬ 
lar hydrogen column density and 12 CO J = 1 — 0 velocity- 
integrated line intensity, and then convert to the units of 
qco via the conversion factor aco (Mg (Kkms^ 1 pc 2 ) -1 ) = 
(2.17 x 10~ 2 °) x Tco(cm -2 (Kkms -1 ) -1 ), as given in 
ISandstrom et all (120131 . in which the ~ 36% contribution 
from helium is accounted for). 

The properties of the minimum-x 2 fits to the 
data for the four clumps are summarized in Table 
[2] and span the range log [n(H 2 )/cm -3 ] = 3.6-3.9, 

log[(Nco/Av)/cm _2 km _1 s] = 15.7-16.5. Three of the 
clumps have minimum-x 2 temperatures 7k ~ 70 K, while 
the fourth appears to be significantly colder: 7k ~ 30 K. 
If the temperatures of the gas and dust are coupled, then 
this range of gas temperat ures is in agreemen t with the 
dual-component dust model dlvison et al.ll2010bh which fits 


the far-IR SED of SMM J21352. We briefly note that our 
measured H 2 densities (n(H 2 ) ) are c oncordant with those 
measured by I Danielson et al.l (l2013l l in their analysis of 
optically-thin 13 CO and C 18 0 emission, while our tempera¬ 
tures are considerably lower than their range of best-fitting 
temperatures (7*, = 140-200 k). However the uncertainties 
on these measurements are large, and at least for clumps 
X, Z\ and Z 2 , our measu red temperature s are w ithin the 
Ax 2 = lcr range quoted bv lDanielson et all (l2013lb Only for 
component Y is our best-fitting temperature formally out¬ 
side their range of allowed temperatures (Tk = 110-200 k). 
The most likely explanation for this apparent disagreement 
is that although the spatially-identified 12 CO J = 1 — 0 gas 
clumps have similar velocity offsets (relative to the sys- 
temic velocity o f SMM J21352) to those clumps studied in 
iDanielson et al.l J2013lb the differing methods by which the 
two sets of clumps have been identified mean that we may 
not be looking at exactly the same material in both cases. 
The low excitation requirements of the 12 CO J = 1 — 0 line 
make it ubiquitous throughout the cool ISM - thus regions 
with high J = l —0 surface brightness correspond to regions 
of dense, cold gas - while both C 18 0 and 13 CO preferen- 
tially trace regions o f recent, intense star-formation (e.g. 
iMeier fc TurneJ 120041 '). in which we may expect to see an 
elevated gas temperature. 


Since radex computes molecular line intensities for a 
given average volume density, n(H 2 ), we must integrate the 
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Figure 5. 12 CO SLEDs of four marginally resolved clumps in SMMJ21352, with the model SLEDs from radex which best fit the 
data over-plotted as a solid black line, radex models are generated across the parameter space 7k=15-90K (in 16 steps of 5 k), 
n(H 2 )= 10 2—6 cm -3 (in 200 equidistant steps) and Nqo — 4 x 10 18—24 cm -2 (in 200 equidistant steps). The key physical parame¬ 
ters of the best fitting (minimum-x 2 ) radex models are given in Table[2] we recover one clump with gas ki netic temperature = 30 k 
and three clumps wi th = 70 k, in good agreement with the distribution of dust tem peratures measured by lSwinbank et al.l (l2010h and 
llvison et al.l d2010tj) . The source-averaged 12 CO SLED from Daniels on et al.l <1201 if) is shown as a dotted grey line for reference. Also 
shown (as shaded regions) are the ±2<r range of models frornXNaravana n Sz Krumh olz (2014, see also mm which predict the 12 CO 
SLED, given estimat es of Ssfr. — we see excellent a greement between the data and both models for clumps X, Z\ and Z 2 , however note 
that in clump Y the iNaravanan &; Krumhold 120 1 4 ) model over-predicts the measured 12 CO J = 6—5 flux. 



Table 2. SMM J21352 clump properties 


Clump 

T k 

(K) 

log [n(H 2 ) 
/cm -3 ] 

“CO 

(Kkms -1 pc 2 ) 

Afgas 

(xlO 9 M 0 ) 

SsFR 

Mq yr -1 kpc -2 

^radio 

9IR 

A' 

70 

3.9 ±0.4 

< 7.8 

< 5.4 

490 ± 220 

- 

2.5 ±0.2 

Y 

70 

3.6 ±0.4 

< 0.6 

< 6.1 

1650 ± 310 

-0.88 ±0.07 

2.3 ±0.1 

Zi 

30 

3.7 ± 0.6 

< 9.4 

< 7.8 

1360 ± 280 

-0.82 ±0.07 

2.2 ±0.1 

z 2 

70 

3.9 ±0.1 

< 16.1 

< 10.6 

510± 220 

- 

2.6 ±0.2 


mean density along a line of sight, f ds, through each cloud 
in order to calculate the column density: Nh 2 = / n(H 2 )ds. 

This conversion from volume density to column density 
of H 2 implies that the best aco is strictly an upper limit; the 
only estimate available for the sizes of the clumps (ds) is the 
constraint that they be no larger than the area covered by 
VLA synthesized beam (~ 100 pc). If the clumps are in fact 
significantly more compact than this, then aco will scale 
downward correspondingly. The resulting limits on aco are 
reported in Table 2. 

The measured upper-limits in aco for the clumps 


are consistent with the range of Ultra-Luminous Infrared 
Galaxy (ULIRG)-like to Milky Way GMC-like values, and, 
in conjunction with estimates on ^CO allow upper-limits 
to the gas masses to be determined (Tablc[2|. Depending 
on whether the total molecular gas mass of SMM J21352 is 
closer to the firm lower or upper limits determined in £13.21 
these clumps may between them contain 10-60% of the total 
molecular gas budget of the system. 
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4-2.1 Comparison with recent theoretical work 

Recently, iNaravanan fc Krumholz ([2014) presented a theo¬ 
retical framework for the excitation of 12 CO in star-forming 
galaxies which draws together numerical simulations of disc 
and merging galaxies with molecular line radiative trans¬ 
fer calculations. Within this framework, it is predicted that 
although in general the shape of the 12 CO SLED is deter¬ 
mined by several physical properties which are observation- 
ally challenging to determine (such as the gas density, tem¬ 
perature and optical depth in the molecular ISM), the fact 
that each of these properties correlates in theoretical models 
in some predictable way with the SFR surface density, Esfr, 
allows the 12 CO SLED to be predicted given a measure of 
Esfr. 

With our latest SMA 870 pm map, we estimate Esfr 
for each of the four clumps. We begin by measuring the flux 
density (Ss 70 /jm) of each clump, and fitting a scal ed versi on 
of the SMMJ21352 far-IR template (Ilvison et al.1 fioiQbl h 
which we integrate between 8-1000 (i m to get Lir. We adopt 
the Lir to SFR conversion factor of iKennicuttl (ll998lb and 
estimate clump sizes based on the angular size of the SMA 
beam (corrected for amplification). 

We use an idl script to generate for each clump 10 5 
Normally-distributed values of Esfr, model, whose central 
value and standard deviation are given by the measured 
Esfr f or t hat clump , an d its_ uncertainty. We then call the 
INaravanan fe Krumholzi (12014) idl code for each Esfr, model 
in the distribution, giving a library of 100,000 12 CO SLEDs 
for each clump. We show the range of model SLEDs whose 
input Esfr, model is within ±2cr of its measured value for 
each clump as a greyscale in Figure[5l along with the ob¬ 
served SLEDs from our deblending work, and the best¬ 
fitting model SLEDs from our radex radiative transfer cal¬ 
culations. We see excellent agreement between both sets of 
mo dels and the data for clu mps A', Z\ and Z 2 , however 
the INaravanan fe Krumholj (j20l4 'l SLED for clump Y di¬ 
verges slightly at J > 4. Hence, to the extent that our low- 
resolution mid- J 12 CO d ata can constrai n the clump SLEDs, 
we End that the INaravanan fe Krumholzi (l2014|j models ad¬ 
equately predict the 12 CO SLEDs of star-forming regions for 
a given Esfr around 75% of the time. 


4.3 The FIRRC and radio spectral index 

It has long been established that in local star-forming galax¬ 
ies, there exists a correlation between the relat ive strengths 
of far-infrared and radio emission (e.g. iHelou et al.1 Il985l : 
ICondon et al.lll99ll ) on scales from CMC complexes up to 
entire galaxies. Furthermore, this correlation is observed to 
hold ubiquitously over several orders of magnitude in lumi¬ 
nosity, gas surface density and photon energy density. 

This FIRRC is most readily unde rstoo d within 
th e framewor k o f ca lorimetry (e.g. iVoelk 1 19891: 
lLacki fc Thompson! hold ). in which UV and optical 
photons produced by young, massive stars are captured by 
dust and re-radiated at infrared wavelengths at a rate that 
is closely matched by that at which the same massive stars 
produce synchrotron emission via supernova shocks when 
they die. Provided the duration of the star formation event 
is greater than the ~ lOMyr time-scale on which the Erst 


supernovae are produced, then the two processes can reach 

equilibrium, giving rise to the observed correla tion. _ 

_The FIRRC is commonly expressed (as in Ilvison et al.l 

boiOai ) via the single logarithmic parameter giR, given as 


9iR = logio 


Sir 


Wm" 2 Hz- 


3.75 x 10 12 W m“ 


Si 


4 GHz 


(3) 


where Sir is the far-infrared (8-1000 pm) flux, and Si. 4 GHz 
is the monochromatic 1.4 GHz radio flux, fc-corrected to the 
rest-frame. 

We investigate this interplay between thermal dust 
emission and non-thermal synchrotron emission in 
SMMJ21352 by making a tapered C-band radio map 
whose resolution is close to that of the SMA map, and then 
convolving this image with a Gaussian kernel to exactly 
match the SMA resolution. Next, we regrid both images, 
binning by a factor 3x in both dimensions, and fc-correct 
the total measured C-band radio flux density in the regrid- 
ded map (using the average spectral index, a ra dio = —0.96) 
to iz r est = 1.4 GHz. Using Lir = (2.3 ± 0.2) x 10 12 Lq as in 
§ 14.2.11 we measure the galaxy-averaged FIRRC parameter 
gin = 2.39 ±0.17. 

Next, we capitalize on the high spatial resolution of our 
new dataset by measuring giR on a (rebinned) pixel-by-pixel 
basis throughout the molecular gas reservoir. We use the re¬ 
binned radio and spectral index maps to create a fc-corrected 
rest-frame 1.4 GHz flux density map, and create a map of 
5*ir by scaling th e observed infrared SED for SMMJ21352 
dSwinbank et al.ll2010l : I viso n et al. I l2010b|) to fit the mea¬ 
sured 870 p -m flux density at each position. We integrate at 
each pixel between 8-1000 fim to obtain a spatially-resolved 
Lir map. In so doing, we explicitly assume that the far-IR 
SED of each resolution element is representable by the same 
mixture of 30 k and 70 k dust as explains the SMM J21352 
SED. We restrict our analysis to pixels which sit within the 
3er contour of 12 CO J = 1 — 0 emission, and for which the 
uncertainty in the spectral index da Ta di 0 ^ 0.2. The median 
value of giR = 2.33 ±0.02 is consistent with that derived for 
SMM J21352 using the total flux densities. 

We present the pixel-by-pixel distribution of a ra dio and 
giR within the molecular gas reservoir in Figure[6] Data are 
shown as coloured dots (colour-coded by 12 CO J = 1 — 0 
flux density), with grey symbols representing measurements 
and limits on Q ra dio and giR in individual SMGs from 
our programme of ALMA observations in the Ext ended 
Chandra Deep Field South (ALESS: iHodee et al.1120131 : 
iThomson et al.ll20l4). Also plotted in Figurc[6] are model 
tracks from Bressan et alJ ((2002), along which starburst 
galaxies are expected to evolve from young to old over a 
timescale of ~ 400 Myr. The pixels at the positions of clumps 
Y and Z\ are highlighted with error bars (which are repre¬ 
sentative of the error bars for all pixels). Clumps X and 
Z 2 lie in the masked region of the spectral index image, for 
which da ra dio > 0.2, and thus we are unable to report their 
«ra,dio. We see that our rebinned pixels in SMM J21352 span 
a narrow range in giR (possibly a direct consequence of the 
assumption of a common dust SED for all pixels), but that 
they exhibit considerable scatter in a ra dio, i ntersecti ng the 
model tracks for evolving starburst galaxies of lBressan et al.1 

(|2002h . 


© 2014 RAS, MNRAS 000.IT1IT51 






























































The ISM on ~ 100 pc scales in a z ~ 2.3 galaxy 11 


CO(1—0) flux (mJy/beam) 


9.9 10.0 10.1 10.2 10.3 

!°g( l irAs) 


1.5 2.0 2.5 3.0 3.5 21.0 21.1 21.2 21.3 21.4 21.5 21.6 21.7 

9ir Iog(Li.4. ghz/WHi -1 ) 

Figure 6. Left: Plot of spatially resolved (i.e. pixel-by-pixel) variations in the radio spectral index, a ra dioi versus FIRRC parameter, q\ r 
(with the error bars shown for clumps Y and Z\ being representative for the majority of pixels). Data points are colour-coded by the 
12 CO J — l — 0 flux density from the smoothed VLA map. We see a tendency for dense gas (high J = l — 0 flux; lighter-coloured points) 
to be co-spatial with regions of higher q\ r (i.e. enhanced infrared emission relative to radio), which is expected if the densest gas is 
co-spatial with the youngest star-forming regions which have not yet produced significant non-thermal radio emission from supernovae. 
For comparison, we plot (in grey) the results of a similar analysis of individual ALMA-observed SMGs from t he on-going ALESS 
programme (data taken from lThomson et al.ii2014l l. Also plotted are model tracks for evolving starburst galaxies llBressan et al JI2002L 
see § !5.2l for details), with red boxes denoting the predicted ages of the galaxies/material they contain. Top-Right: Spatially-resolved q \r 
versus log(L|u) and Bottom-Right: q\ r versus log(L ra d; 0 ), from which we see that that spatial variations in q\ r across the surface of 
SMM J21352 are driven predominantly by the radio morphology of the source. 
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4.4 The Schmidt-Kennicutt Law 

A critical ingredient in any attempt to understand galaxy 
formation and evolution is a proper characterization of the 
processes by which diffuse baryonic matter comes to be 
locked up in stars. Unfortunately, the spatial scales on which 
star formation actually occurs remain beyond observational 
reach for all but the most nearby galaxies, and thus efforts to 
characterize these processes in distant galaxies are limited to 
the application of simplified recipes, or star-formation laws. 

It has long been claimed (e.g. iGreve et al . IHqqI) that 
SMGs follow a Schmidt-type star-formation law, Esfr oc 
SgL (where Esfr and E gas are the SFR and molecular 
gas surface densities, respectively, and N is the power- 
law index), albeit possibly with an enhanced SFR per unit 
molecular gas relative to the population of normal star¬ 
forming galaxies; this has been inter preted as ev idence 
of a dual star-formation l a w at work llDaddi et al.l 120081 . 
l2010bl : iDannerbauer et al.l 120091 ). Und er this paradigm , 
photometrically-selected BzK galaxies (IDaddi et al.l 120041 ) 
typify the main sequence of ongoing star-formation activ¬ 
ity, fuelled by cold mode accretion from the surrounding 
intergalactic medium, while SMGs form stars in shorter- 


lived, hy per-effic ie nt sp urts due to major merger activity 
fe.g. lGenzel et al.l20ld ). However, many of the studies which 
have been claimed to corroborate this model have relied on 
observations of the mid-J 12 CO lines, which, while brighter 
(and hence easier to observe) offer a view of the molecular 
gas reservoir that is biased by their excitation requirements 
(i.e. denser, warmer gas, which may not be representative of 
the bulk of gas reservoir). Such observations are reliant on 
an assumed gas excitation model, fr om which aco (and ul¬ 
timately M gas ) is inferred (see llvison et al . I l201li . for further 
discussion of this point). 

The high angular resolution of our new VLA maps 
for SMMJ21352 allow us to evaluate the spatially-resolved 
Schmidt-Kennicutt relation in a high-redshift star-forming 
galaxy in 12 CO J = l — 0, without the need to assume a gas 
excitation model. Tapering the 12 CO J = l — 0 uv data and 
convolving the tapered image with a Gaussian kernel, we 
create a map of molecular gas emission with the same res¬ 
olution as the updated 870 fim map, and measure Lco(i-o) 
in each of the dense gas clumps. 

We use estimates of the SFR of each clump from 
§ 14.31 and convert L' co ^ 1 _ 0) to M gas via a suitable choice 
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of aco- The lowest qco produced by our radiative trans¬ 
fer modelling of the individual clumps (Tabled is aco = 
0.6 M© (Kkins^ 1 pc 2 ) -1 , in clump Y. For A', Z\ and Z 2 , the 
qco upper limits are consistent with the galaxy-averaged 
Qco ^ 3.5 ( S 13.2I) . and so we adopt this limit for these three 
clumps. It will not be possible to significantly improve upon 
this assumption until even higher resolution 12 CO observa¬ 
tions allow us to better estimate the sizes of the clumps, 
and thus better constrain the conversion from H 2 column 
density to number density. 

Acknowledging the limitations imposed by these various 
assumptions, the measured E gas and Esfr of the clumps are 
consistent with an extrapolation of the Schmidt-Kennicutt 
Law (Figure[7]), making them among the densest star¬ 
forming environments seen at high-redshift. SMMJ21352, 
too, unsurprisingly appears to obey the Schmidt-Kennicutt 
Law, albeit at lower surface density. The positioning of the 
clumps on the EsFR/E gas plot sees them slightly offset from 
the star-formation law in “normal” galaxies, but is subject 
to the uncertainties in the values of aco used to derive them. 


5 DISCUSSION AND CONCLUSIONS 

5.1 New constraints on the spatially-resolved 
12 CO SLED 

We have presented new, high resolution VLA observations 
of 12 CO J = 1 — 0 and C-band 4-8 GHz continuum emis¬ 
sion in the strongly lensed but otherwise typical SMG, 
SMM J21352—0102. There being no evidence of flux having 
been significantly resolved out on long baselines, we con¬ 
sider the new 12 CO J = 1 — 0 flux density, luminosity and 
gas mass estima t es - w ithin 5% of the values measured by 
[Danielson et al.l d 201 it ) based on single-dish GBT measure¬ 
ments - to be minor revisions to those figures. 

Taking full advantage of the ~ 0.3” resolving power of 
the VLA’s BnA configuration in A'a-band, we have identi¬ 
fied four dense gas clumps in the J = 1 — 0 data cube, with 
counterparts mirrored about the lensing critical curve. The 
velocity offsets of the 12 CO J = 1 — 0 line profiles of these 
clumps relative to the systematic velocity strongly suggest 
that the y are related to the clumps A , Y, Z\ and Z 2 iden¬ 
tified by iDanielson et al.l (1201 ll . 12013 ), whose analysis was 
based on the kinematic decomposition of the 12 CO spectra 
in to multiple Gaussian components. 

We have investigated the 12 CO excitation conditions 
within these ~ 100 pc sub-structures, determining the likely 
gas temperatures and molecular hydrogen volume densities 
through modelling of their 12 CO SLEDs. This has allowed us 
to place plausible upper-limits on the value of the 12 CO-to- 
H 2 conversion factor aco for each of the clumps, which fall 
between those normally taken for ULIRGs and Milky Way 
GMCs. Depending on whether we adopt the firm upper or 
lower limits on the total molecular gas in SMM J21352, we 
find that these clumps may contain between 10-60% of the 
total molecular gas in the galaxy. 




\JM qP c- 2 ) 

Figure 7. The Schmidt-Kennicutt plot, relating star formation 
and the availability of molecular gas. Local discs, LIRGs, ULIRGs 
and high-red s hift s tar-forming ( BzK ) galaxies are taken from 
iGenzel et al.l (120101) . SMGs with recent 12 CO J = 1 — 0 mea¬ 
surements, including the GN 20 object of lCarilli et al.l feOlOl) are 
shown. Top: The observed S-K relation. After correction for lens¬ 
ing, SMM J21352 is seen to be an intrinsically typical SMG, while 
the analysis of the clumps extends the plot at high redshift to lu¬ 
minosities hitherto well below the confusion limit. Bottom: The 
surface density S-K relation, in which we convert the observed 
12 CO J — 1 — 0 luminosities in to gas masses using a standard 
ULIRG-like aco = 0-8 for the whole galaxy (to facilitate com¬ 
parison with the literature) and qco — 0.67 for the clumps, based 
on the radiative transfer calculations presented in $ 14.21 On both 
plots, the error bars shown are calculated from the signal-to-noise 
of the SMA 870 pm and 12 CO J = 1 — 0 images, on the assump¬ 
tion that the clump sizes are well-known. If the clumps are in fact 
much smaller than the VLA beam, then both Esfr and E gas are 
raised proportionally, but with an additional raising of E gas due 
the lower aco that this increase in surface density would imply. 
The dense clumps in SMM J21352 populate the high-density tail 
of the so-called “merger sequence”, however systematic uncertain¬ 
ties in aco in both BzK galaxies and SMGs make the distinction 
between these two sequences unclear. 
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5.2 The FIRRC/spectral index relation as an age 
estimator 

Combining our updated SMA 870 p ,m map with brand-new 
4-8 GHz C-band continuum observations of SMMJ21352, 
we measure giR, and hence test whether SMM J21352 obeys 
the far-infrared/radio correlation. We measure giR = 2.39 ± 
0.17 for SMM J21352, with individual pixels spanning a nar¬ 
row range giR ~ 2.1-2.6, which may have some interesting 
implications for th e age di stribu tion of material in the ISM. 
According to the iBressan et alJ 12002^ models of starburst 
galaxies, slight deviations from the FIRRC are expected for 
star-forming systems whose luminosity-weighted age is suf¬ 
ficiently young. Whereas the dust heating which produces 
thermal infrared emission is nearly instantaneous with the 
onset of star-formation, the radio spectrum is at early times 
comprised only of (relatively weak, flat-spectrum) thermal 
Hu emission, and only at later times do the first supernovae 
- which produce (steep spectrum) non-thermal synchrotron 
radio emission - occur. The net effect of these differing time- 
scales is that both a ra dio and giR for a typical galaxy may 
go through three key phases: (i) the proto-starburst phase, 
during which dust-heating by young, massive stars rapidly 
elevates giR, and the radio brightness is low, dominated by 
thermal emission; (ii) the evolving starburst phase, during 
which time the heating of dust by low-intermediate mass 
stars continues to elevate Sir, but is out-paced by the much 
more rapid elevation in L ra dio due to synchrotron emission 
originating in the first supernova remnants, and; (iii) the 
post-starburst/quiescent phase, which is much longer-lived 
than either of the preceding phases, and occurs as the star- 
burst event draws to a halt, and the heating of dust and the 
production of supernovae reach equilibrium, whereupon giR 
and a asymptote. 

If the iBressan et all ( 2002 .1 models accurately describe 
the far-infrared and radio emission in SMM J21352, then it 
is notable that we see so few pixels with giR < 2.2, occupying 
the “evolving starburst” phase. The observed distributions 
in Q ra dio and gm appear to overlap the model tracks at both 
the young, proto-starburst (< 30Myr) and the old, post- 
starburst (> 200 Myr) ends. While it is certainly possible 
that the ISM in SMMJ21352 contains material of different 
ages, we consider the bi-modal age distribution implied by 
the scatter of the data about the model tracks to repre¬ 
sent an implausible scenario. Given the representative error 
bars shown in Figure[6] we cannot rule out the possibility 
that all the pixels have Q ra dio an d gm values that lie at the 
young ends of the iBressan et all d2002i l evolutionary tracks; 
indeed, given th e high instantaneo us SFR in SMM J21352 
(~ 400 Mq vr~ 1 : llvison et al .20103) , we may expect to mea¬ 
sure a low average age for material in the ISM. That the data 
are plausibly consistent with the entire ISM having an age 
< 30 Myr (there bei ng no data covering the “evolving star- 
burst” region of the IBressan et alj|2002l . model tracks with 
giR < 2.2) may indicate that SMM J21352 is a near-pristine, 
proto-starburst galaxy, whose high SFR has been rapidly 
triggered by a recent interaction. 

An alternative explanation for the spread in giR seen 
within SMM J21352 is that it may be due to differing diffu¬ 
sion timescales and mech anisms for cosm i c ray electrons and 
UV photons in the ISM. iMurphv et al.1 (2006;) find similar 
offsets in the logarithmic 24 /im and 70 ^tm-to-radio surface 


brightness ratios for ~ kpc regions within the discs of four 
nearby star-forming galaxies. These offsets disappear when 
the infrared maps are smoothed with exponential kernels, 
which attempt to repli cate empirical “leaky box” cosmic 
ray electron models (e.g. lBicav fc Heloull 1990l : [Murphy et ali 
1200614 . In these models, cosmic ray electrons escape from 
the young proto-stellar environment more quickly than do 
the thermal infrared photons, which can remain bound up 
in dust knots for > 100 Myr before eventually escaping via 
a random walk process. This leads to the radio emission 
having a correspondingly larger scale length (~ 1-5 kpc, cf. 
~ 100 pc for the far-IR) and less clumpy structure than the 
infrared emission. Similar kpc-scale offse ts in the FI RR C are 
seen in clumps in nearby spiral galaxies dDumas et al. liml; 
iBasu etld1l2012l 'l. 

There are thus two competing effects which may ex¬ 
plain the pixel-to-pixel distribution in giR and a ra dio within 
SMM J21352, which operate on different time-scales and im¬ 
ply opposing estimates for the ages of the ISM: (i) that the 
offsets are due to the ~10 Myr lag between dust heating and 
the onset of the first supernovae, and reveal the existence of 
a young, near-pristine ISM; or (ii) that the offsets are due 
to the cosmic ray electron diffusion time-scale being signifi¬ 
cantly shorter than that of the infrared photons, which can 
remain bound up in dust knots for > 100 Myr before eventu¬ 
ally escaping via a random walk process, implying an evolved 
ISM. The present data do not have the ability to fully distin¬ 
guish between these two scenarios, however we can hope to 
revisit this dichotomy in future if we can obtain an estimate 
of the time-scale on which radiation pressure and mechani¬ 
cal feedback could conspire to disrupt the dense gas clumps 
in SMMJ21352. Attempts to understand of the details of 
how galactic velocity fields are coupled to the turbulent dy¬ 
namics of their constituent clouds - a necessary ingredient of 
such an age estimate - do, however, rem ain an area of ongo¬ 
ing r esearch even in the local Universe dKrumholz fe Dekell 
l2Q10l. and references therein). 

From these arguments, we see that the observed varia¬ 
tions in gm for the star-forming clumps within SMM J21352 
are not unexpected, however a more rigorous investigation 
of the spatially-resolved FIRRC will require forthcoming, 
high-resolution ALMA observations, which will allow us to 
create spatially-resolved far-infrared SEDs for each of the 
clumps - and so obtain robust estimates of Lir for regions 
within SMMJ21352. This will eventually free us from the 
assumption that the SMM J21352 far-IR SED can simply 
be scaled to match the 870 /tm flux density of each clump. 

5.3 The star-formation law in 2 ~ 2.3 GMC 
complexes 

We have used the new 12 CO J = 1 — 0 VLA map in con¬ 
junction with the SMA 870 /im map and existing mid-J 
12 CO maps from PdBI to investigate the properties of the 
cold molecular gas inside each clump, and ultimately place 
both the galaxy as a whole, and the clumps on the Schmidt- 
Kennicutt plot. The implied upper-limits on the 12 CO-to-H 2 
conversion factor, aco, derived via radiative transfer cal¬ 
culations are consistent with those previously derived for 
SMGs, though these limits are sensitive to the sizes esti¬ 
mated for the clumps. Within the limitations of these as¬ 
sumptions, we see that SMMJ21352 and the clumps lie 
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within the scatter of the Lm/L' co plot, at the moder¬ 
ate and low luminosity ranges, respectively, in good agree¬ 
ment with the findings an d implications of recent 12 CO 
J = 1 — 0-based surveys (e.g . lHanis_et_^l Il2010l: llvison et ahl 


J = 1 — U- b ased surveys (e.g. Harris et a l.HZUIUI; llvison et alJ 
l201ll : iRiechers et all l201ll : Thomson et al] 120121 ). While 


these clumps appear to be marginally offset from the “nor¬ 
mal” star formation law claimed for BzK galaxies on the 
SsFn/Sgaa surface density plot, we note that: (i) M gas is 
usually very poorly constrained in BzK galaxies due to the 
reliance on an assumed brightness temperature ratio to pre¬ 
dict the 12 CO J = 1 — 0 line flux from mid-J observations, 
and on the su bsequent app lic ation of a standard aco to all 
BzK galaxies llDaddi et al.ll2010ah . leaving scope for consid¬ 
erable variation in the values of Esfr, an d E gas in individual 
BzK galaxies that is not reflected in the scatter in Figure]?] 
(ii) using a lower aco for SMM J21352 - which is allowed 
by our radiative transfer modelling - would imply a higher 
gas mass surface density (at a fixed Esfr), and so would 
shift the clumps to the right in Figure]?] and in to better 
agreement with a single star formation law. 


5.4 The relationship between 12 CO excitation and 
star formation 

Finally, we have used our suite of multi-wavelength data 
to test the predictio n of a recent theoretical model 
dNaravanan fe Krumhola feoi-lJ) which claims that although 
the 12 CO excitation conditions in galaxies are determined by 
a morass of difficult-to-observe quantities (including tem¬ 
perature, density and optical depth), each of these corre¬ 
lates in some predictable way with the SFR surface density, 
Esfr, such that a measurement of Esfr can yield a predic¬ 
tion for the shape of the 12 CO SLED. Using SMA 870 Mm 
observations of SMM J21352, in conjunction with the well- 
sampled far-IR SED of the galaxy, we measure Lir, SFR 
and Esfr for each of the dense molecular gas clumps iden¬ 
tified in 12 CO J = 1 — 0, and create a prediction for the 
12 CO SLED. We find that in three of the four clumps, the 
predicted SLEDs are in good agreement with those derived 
from multi-J 12 CO imaging, providing tentative support for 
the iNaravanan fe Krumhold (|2014| j model of 12 CO excita¬ 
tion in star-forming galaxies. 
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